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ABSTRACT 

Context. Fourier transform (or lag) correlators in radio interferometers can serve as an efficient means of synthesising spectral chan- 
nels. However aliasing corrupts the edge channels so they usually have to be excluded from the data set. In systems with around 
10 channels, the loss in sensitivity can be significant. In addition, the low level of residual aliasing in the remaining channels may 
cause systematic errors. Moreover, delay errors have been widely reported in implementations of broadband analogue correlators and 
simulations have shown that delay errors exasperate the effects of aliasing. 

Aims. We describe a software-based approach that suppresses aliasing by oversampling the cross-correlation function. This method 
can be applied to interferometers with individually-tracking antennas equipped with a discrete path compensator system. It is based 
on the well-known property of interferometers where the drift scan response is the Fourier transform of the source's band-limited 
spectrum. 

Methods. In this paper, we simulate a single baseline interferometer, both for a real and a complex correlator. Fringe-rotation usually 
compensates for the phase of the fringes to bring the phase centre in line with the tracking centre. Instead, a modified fringe-rotation 
is applied. This enables an oversampled cross-correlation function to be reconstructed by gathering successive time samples. 
Results. Simulations show that the oversampling method can synthesise the cross-power spectrum while avoiding aliasing and works 
robustly in the presence of noise. An important side benefit is that it naturally accounts for delay errors in the correlator and the 
resulting spectral channels are regularly gridded 

Key words, instrumentation: interferometers - techniques: interferometric 



1. Introduction 

The observing band of radio interferometers frequently needs to 
be split into sub-bands, either for spectral line observations or 
to reduce the effects of chromatic aberration. Fourier transform 
correlators offer an efficient method for dividing the observa- 
tion band. In this scheme, the signal from the two arms of the 
interferometer are correlated at discrete delay steps, making di- 
rect measurements of the cross-correlation function (see Fig.!!). 
Taking the Fourier transform of the cross-correlation function 
gives the complex cross-power spectrum. For a signal of band- 
width Av, Nyquist sampling theorem requires the signal to be 
sampled at time intervals of l/(2Av) to avoid aliasing. But the 
cross-correlation function of a band-limited signal extends over 
an infinite delay range. So the Nyquist sampling theorem holds 
true only if we sample the cross-correlation function over an 
infinite delay range. Clearly this is not practical so the cross- 
correlation function is sampled over a finite range. This results 
in a recovered signal spectrum with tapered band edges which 
will overlap with its images in the spectral domain. This over- 
lap causes aliasing and will corrupt the recovered cross-power 
spectrum. 

As an illustration, the plots in Fig. |2] show snapshots of the 
cross-correlation functions and the recovered spectra. The left- 
most plots (case 1) show the cross-correlation function when the 
source transits an east-west baseline (that is when the path dif- 
ference between the two arms of the interferometer is zero). The 
circles represent the measurements at discrete delay steps of the 
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correlator. The next plots to the right (case 2) show the cross- 
correlation function a short while later. The right-most plots 
show the spectrum calculated by taking the Fourier transform 
of the discrete measurements. The filled grey circles indicate the 
spectrum at transit (case 1) and the dark circles are for case 2. In 
the top row of Fig. [2^, the cross-correlation function was sam- 
pled at 16 delays. This critically samples the Av = 6 GHz signal 
bandwidth but undersamples the 0-12 GHz basebandsignal, so 
the positive and negative halves of the 6-12 GHz bands lie side- 
by-side. In case 1, the cross-correlation measurements trace out 
a delta function and give the characteristic flat spectrum. But 
in case 2, the amplitudes are perturbed. Ideally, the amplitudes 
should not change so this points to a fundamental problem. This 
could be worked round by oversampling the baseband signal at 
64 or more delays as illustrated by Fig. |2j3. Now the two sets 
of amplitudes at different times match perfectly. Oversampling 
introduces a buffer between the two halves of the passband and 
also between their spectral images. This suppresses aliasing by 
reducing the overlaps between the signal and image bands. But 
clearly, sampling at 64 delays in hardware is not practical. 

This is manifested in the Fourier Transformed data (the spec- 
trum) as temporal modulations in both the amplitude and phase 
(see Fig. [3]l. In aliased signals, the noise components between 
the channels will also be correlated so the individual channels 
cannot be strictly treated as independent measurements. These 
effects are particularly pronounced in the edge frequency chan- 
nels as seen in Fig.|3]so these channels are usually rejected. This 
may be acceptable in systems with tens to hundreds of channels. 
But some correlators, particularly broadband analogue Fourier 
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Fig. 1. a) Top plot: The cross-correlation function measured by 
a real correlator (black) and the envelope (grey). This is a graph- 
ical representation of the cross-correlation function in Eq. ([T) 
with the geometrical delay Tg = 0. The open circles indicate the 
position of the detectors in instrument delay t-,. Bottom plots: 
The cross-power spectrum recovered by applying the Fourier 
transform on the cross-correlation function. We have assumed 
the telescope parameters in Table [1] b) The in-phase (real) and 
quadrature (imaginary) cross-correlation functions sampled by a 
complex correlator. The cross-power spectrum recovered in this 
case is single-sided. 



transform correlators with channels of order 10 (for example, 
[Li et al. 2004; Roberts et al. 2007; Holler et al., i2007l hereafter 
H07), the loss constitutes a significant portion of the total band- 
width. 

In principle, there are three ways to suppress aliasing: 

1 . Increasing the delay range over which the cross-correlation 
function is sampled. The spectral channel width will be nar- 
rower and the edges of the passband will be sharper. So the 
overlap between the signal spectrum and its images will be 
narrower. The edge channels will still have to be rejected but 
it will be a smaller portion of the whole data. 



2. Oversampling the cross-correlation function at finer delay 
steps as already illustrated with Fig.|2] 

3. Reducing the bandwidth so that it is not critically-sampled. 

Implementing either of the first two modifications in hardware 
would be costly and the third approach would lose sensitivity. 
In Sect. [3] we propose a software-based approach to oversample 
the cross-correlation function and avoid aliasing. This is based 
on the well-known notion that the spectrum of a source can 
be obtained from the cross-correlation function measured by a 
drift scan. We can reconstruct the cross-correlation function us- 
ing a combination of source tracking, path compensation and 
fringe-rotation. In Sect. IH we will illustrate an application of 
this technique with simulations. The system requirements are; 
(1) Individually-tracking antennas and (2) discrete-delay path 
compensation. We have modelled an analogue correlator here 
but the principles could equally apply to digital correlators. 

A number of groups recently reported broadband analogue 
Fourier t ransform correlators sufferi ng fr om delay errors (fo r 
example. [Harris & Zmuidzinasll2001l: H07 ; [Roberts et al.ll2007h . 
Simulations by H07 showed that delay errors make the effects 
of aliasing worse. Although t he delay errors can be calib rated a 
sample at a time (recently bv [Harris & Zmuidzinasll200ll ; H07), 
the method described here will be a natural way of accounting 
for delay errors. This side benefit is perhaps as significant as the 
alias suppression aspect of the method. The resulting spectral 
channels are also regularly gridded at the desired frequencies. 

However there are a number of issues that must be consid- 
ered in a practical system and these are discussed in Sect. 14.31 
This paper is a follow up to H07 where we described the de- 
velopment of a broadband 6-12 GHz Fourier transf orm cor- 
relator for the Arcminu te Microkelvin Imager (AMI; [Kanekol 
[2006[ ; [Zwart et al.[[2008[) . AMI is a new interferometer designed 
to survey for clus ters of galaxies by exploitin g the Sunyaev- 
Zel'dovich effect JSunvaev & Zerdovich[[l972[) . The effects of 
aliasing in the correlator used in this instrument are discussed 
further in H07. Before going into a detailed discussion of the 
oversampling method, we will first give a brief overview of how 
the data is processed in a conventional interferometric system. 



2. Conventional interferometry 

As the antennas of an interferometer track a source, the geo- 
metric path difference (rg) between the two arms needs to be 
corrected (Fig.H]). One approach called path compensation is to 
insert discrete lengths of delay lines (rpc). Usually, the signal 
at the observation frequency Vrf is downconverted by mixing it 
with the local osciallator (LO) at Vlq. This converts the signal to 
a lower intermediate frequency (IF) at Vip. The path compensa- 
tion is usually inserted in the IF and we will assume that this is 
the case. It is possible to eliminate the path compensation sys- 
tems all together by co-mounting the antennas. But our spectral 
synthesis method cannot be applied to such a system and must 
have individually-tracking antennas. Individually-tracking aten- 
nas also give better control on systematic errors because non- 
astronomical noise and other contaminating signals can be re- 
moved by fringe-filtering. 

An interferometer is a device that measures the cross- 
correlation function of the signals received by a pair of antennas. 
For a source at the tracking centre, this is given by 



sin;rAv(Tg-Ti) 
Ri = Rq cos 

7rAv(Tg - Ti) 



[27r(vLorg + v,F(Tg - Ti) -H (^Lo)]- (1) 
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Fig. 3. Simulated time-stream data of the amplitude and phase 
for each spectral channel. The cycles are caused by aliasing and 
are worst in the edge frequency channels. This simulation is for 
a real correlator with parameters outlined in Table [T] Similar 
effects are also seen in the complex correlator. 




basis of Fourier transform correlators that we will be discussing 
in this paper 

For brevity, we will define the residual delay at each correla- 
tor delay as 



Ar = Tg - Ti = Tg - (Tpc + rig). 



(3) 



We can then re-express Eq. ([TJ using At and this will be the 
standard form we will use in the rest of the paper; 

sin.7rAvAT r / \i 

Ri = Ro — -r-j — cos[27r(vLorg ± v.fAt + 0^0 jj- (4) 



2. 1 . Real and complex correlators 

Broadly speaking, there are two architectures for Fourier trans- 
form correlators. We will call these the real and complex corre- 
lators. The spectral information they give are the same but the 
way they sample the cross-correlation function is different. A 
more detailed discussion on Fourier transform correlator designs 
can be found, for example, in H07. The real correlator samples 
at 2N delay steps and will produce A^ sub-bands. The complex 
correlator makes two measurements at each of the A^ delay steps; 
the in-phase component (like in the real correlator) and also the 
quadrature component. The quadrature component is detected 
by inserting a -90° phase shift into one arm. For the quadrature 
component, the cosine in the cross-correlation function (Eq. (|4|) 
is replaced with a sine; 



sin .TT AvAt 

^2 =^0 7— A — sin 

ttAvAt 



[27r(vLoTg + VifAt + 0Lo)]. 



(5) 



We have also substituted in the residual delay At defined in 
Eq. (O for clarity. The complex correlator samples the in-phase 
and quadrature components at A^ delay steps. The two compo- 
nents can be conveniently expressed by a sum of the two orthog- 
onal signals, 



R.^Ri+ jRi 



(6) 



Fig. 4. A schematic o f an interferometer system based on or written as an exponential; 
[Thompson et al.l(l200lb . The geometric delay (Tg) between a pair 
of antennas is corrected by the path compensator delay (xpc) in 
the IF. The coiTelator samples the cross-correlation function at a 
number of delay steps (rig). 



sin;TAvAT r / \i 

Rx = ^0 — — r — exp[27rj(^VLoTg + v,pAt + 0loJJ- 



(7) 



The coefficient Rq is proportional to the flux density of the point 
source and also includes factors like the gain of the detector cir- 
cuit. The plus sign in the + term applies to upper side-band sys- 
tems and the minus sign is for lower side-band systems. 0lo is 
the phase of the local oscillator and it can be removed by phase 
calibration, at a later stage, t; is the total instrument delay dif- 
ference in the IF inserted by the path compensator (Tpc) and the 
differential delays inserted by the coiTelator lag (rig); 



Ti 



-HTlg. 



(2) 



The cross-correlation function (Eq. ([T])) plotted against the cor- 
relator delays rig (at v,f) consists of a cosine term bound by a 
sine envelope (see Fig.[T^). The Fourier transform of the cross- 
correlation function gives the spectrum. So if we can sample the 
cross-correlation function sufficiently at regular correlator de- 
lays Tig, we can measure the spectrum of the source. This is the /ft = exp -Inuvi^oTg + v,fAt) 



The delay steps of the Fourier transform correlator are deter- 
mined by the total bandwidth of the correlator For the critically- 
sampled real correlator, the delays are in steps of T = l/(2Av). 
The complex coiTelator samples at only half-Nyquist rate (T = 
1/Av) because the cross-power spectrum of a complex cross- 
correlation function is single-sided (see Fig.fTh). 

2.2. Fringe-rotation 

Conventional fringe-rotation transforms the data so that the 
phase centre maps to the tracking centre. This is done by com- 
pensating for the cosine component of Eq. (|4|. Fringe-rotation 
can be applied either in hardware by controlling the phase of 
the local oscillator (0lo) or in software after recording the data. 
We will assume the latter case. The data must be sampled fast 
enough to avoid fringe-washing. For the complex correlator, the 
fringe data (Eq. ^) is multiplied by 



(8) 
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This fringe-rotation factor can be calculated exactly from a com- 
bination of the geometry of the baseline (antenna pair), the posi- 
tion of the source and the frequencies Vlo and v,f. Fringe-rotation 
will stop the phase-wrapping of the the fringes to give 



sin ;tAvAt r n 

Rxfb = ^0 — -T—r — exp[27rj0LoJ- 



(9) 



The Fourier transform of the stopped fringes gives the spectral 
data or visibilities. The visibilities can then be used to map the 
source. 

For the real correlator (Eq. (HJi and as illustrated by Fig [1^), 
we cannot directly rotate the fringes in software because we do 
not have both the amplitude and phase information at each de- 
lay step. So the data must first be Fourier transformed across 
the delays into spectral channels (which is the usual operation 
for Fourier transform correlators). Each spectral channel is then 
fringe-rotated. The data from the complex correlator can also be 
fringe-rotated after the Fourier transform. This may be advan- 
tageous, for example when the fractional bandwidth (Av/vrf) is 
high and the centre frequency used for fringe rotation is less con- 
strained. 

2.3. Spectral synthesis and mapping 

In a conventional signal processing pipeline, the data recorded 
at each time sample are Fourier transformed across the delays to 
give the spectral data at each sub-band. For example, the contin- 
uous Fourier transform of Eq. ^ gives a rectangular spectrum 
as expected. The amplitude is calibrated against a flux calibrator 
source and corrected for changes in the system temperature (in- 
cluding effects like airmass). The phase is also calibrated against 
a phase calibrator source. This will remove phase offsets like <pi^o- 
In the ideal case, an unresolved source at the tracking centre will 
produce visibility data with a constant amplitude and zero phase. 
The successive visibility samples are time-averaged to reduce 
the data load and gridded to the two-dimensional aperture plane. 
This gives the spatial frequency components of the map. The co- 
ordinates of the visibility data in the aperture plane depends on 
the baseline geometry relative to the source and the RF of the 
sub-band. When all the visibilities have been gridded and ap- 
propriately weighted, the aperture plane is Fourier transformed 
to the map plane. The result is the dirty map - the convolution 
of the map with the synthesised beam. The synthesised beam of 
an interferometer usually has large sidelobes so the dirty map 
can be deconvolved using standard techniques such as CLEAN. 
Parts of the steps for reducing the data are summarised in Fig.|5] 



3. Oversampling method 

As we already showed in Sect. [T] it is possible to reduce alias- 
ing by oversampling the cross-correlation function above the the 
Nyquist rate. This inserts a buffer zone in the spectral domain be- 
tween the signal spectrum and the images. Using the new method 
that we will describe here, we could achieve this without increas- 
ing the number of detectors in the correlator 

The source is tracked as usual and the path differences be- 
tween the antennas are compenstated in discrete steps using path 
compensators. Normally, the data are fringe-rotated to align the 
phase centre with the tracking centre. This stops the fringes. 
Instead, we apply a modified fringe-rotation (Sect. 13. Il l and al- 
low the fringes to vary. When the path compensator steps, the 
data since the last step are collected together. The varying fringes 
trace out the cross-correlation function in Eq.|2]with Tg - 0. The 



function is densely sampled so the Fourier transform (Sect. 13.21 ) 
introduces negligible aliasing]. 



3.1. Modified fringe-rotation 

The aim of the method is to build up an oversampled set of 
data points for the cross-correlation function in At. But this is 
complicated because the exponential term in the response of the 
complex correlator (Eq. (|7])) is also dependent on Invi^oTg- It is 
possible to remove just this term from the equation by applying 
a modified fringe-rotation function instead; 

/mfr = exp[-2;rjvLoTg]. (10) 

The resulting fringe-rotated data will be a function of At; 



sin;rAvAT r / \n 

^x/mfi- = ^0 — T—r — exp[27rj^+v,FAT + 0loJJ- 



(11) 



By collecting consecutive time samples, we can build up an over- 
sampled measurement of the cross-correlation function. We as- 
sume that the underlying signal is not changing over the time it 
takes to sample the whole cross-correlation function and we will 
return to this issue in Sect. 13.21 

In contrast, the real correlator is not very well suited to the 
oversampling method. The fringes cannot be directly fringe- 
rotated so it first needs to be Fourier transformecfl The spec- 
tral channels are fringe-rotated by the modified factoiQ/mfi from 
Eq. ( fTOl ) as before and then inverse Fourier transformed back; 



Ri = FT-'[/„rf,FT[7?i(Av,AT,Tg)]]. 



(12) 



Note that this gives the rotated fringes as a function of At. In 
addition to the computational burden, the Fourier transform step 
will cause aliasing. We will investigate this phenomenon through 
simulations in Sect. 14.2.11 This shortcoming for the real corre- 
lators may be overcome with a hardware fringe-rotation system. 
This approach also has the benefit of slowing down the fringes so 
that the data can be sampled at a lower rate. If fringe-rotation is 
to be applied in software, complex correlators are preferable be- 
cause direct fringe-rotation is faster and much cleaner However 
broadband -90° phase shifters can be difficult to design. 

So far we have assumed an ideal square passband. In re- 
ality, the source's spectrum may not be flat and the instru- 
ment may have a sloping passband. These factors may dis- 
place the centre-frequency and degra de the system's sensitiv- 
ity. But following methods outlined bv^Thomp son & D'Addarid 
(1982); Thompson et al. (2001), we calculated that the losses 
from fringe-rotation are negligible compared to losses from the 
sloping spectrum. 



' The oversampling method shares some analogies with drift scans. 
In a drift scan, the source is tracked but the path compensator is kept 
fixed. The spectrum can then be estimated by taking the subsequent 
Fourier transform of the time-series data. 

- In signal processing, window functions are often applied to reduce 
sidelobes. But a window function should not be applied at this stage 
because the spectral resolution needs to be as narrow as possible to 
avoid aliasing - aliasing from the overlapping main lobes is a more 
acute problem. 

^ The conjugate relationship between the positive and negative spec- 
tra in a real correlator means that when applying fringe-rotation to a real 
correlator, the conjugate of Eq. dlOt may have to be applied to one-half 
of the spectrum 
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3.2. Spectral synthesis 

We now have samples of the fringe-rotated cross-correlation 
function against At. We now divide the time-series data between 
path compensation changes so that in each PC data block, the 
cross-correlation function is completely sampled over the re- 
quired range in At. The PC data block is then gridded and each 
pixel is weighted by the system temperature at the corresponding 
sample time. The gridded PC data block is then Fourier trans- 
formecQ to give the oversampled spectrum. The wide spectral 
buffer between the signal and the images reduces aliasing. The 
synthesised spectrum over the passband is extracted, calibrated 
and gridded to the aperture (mv) plane as described before. This 
process is repeated for each PC data block. We have summarised 
the steps as a flowchart in Fig. |5] and we will describe some ex- 
amples through simulations in the next section. 



Data collection: 




Flag slow fringes 
and corrupted data 



Collect data until 

path compensator 

switches 



Oversampling Method: 

} 

This PC data block 
sufficiently spans the CCF 



Standard 
Procedures: 



1 



Apply modified 
fringe-rotation 



data is now 
In AT-space 



Grid data 
(samples the CCF) 



Fast Fourier 
Transform the CCF 



Select spectral data 
over the physical passband 



No 




Calibrate amplitude 
and phase against 
calibration sources 



Grid to aperture 
(uv) plane 



FFT to produce 
a dirty map 



CLEAN 



T 

^Map^ 



Fig. 5. A flowchart summarising the oversampling method. 

The period of the PC data block depends on the baseline and 
will vary over the course of an observation. A fundamental re- 
quirement is that the PC data block period n eeds to be short 
enoug h to avoid time-average smearing (see [Thompson et al.l 
1200 lb . So the periods of the PC data blocks T^^ must satisfy 

Tp, < (£//D)/a, (13) 

where d is the antenna diameter, D is the baseline and Qe is the 
angular speed of the Earth's rotation. The period of a PC data 
block for a complex correlator is given by: 

T^ = (v,p/Ay)/vf, (14) 

where Vf is the fringe frequency. For real correlators, where the 
shortest path compensation bit is shorter, T^. is half the above 



'' A window function could be applied to the gridded data before the 
Fourier transform step but it needs to be chosen with care; we discussed 
in H07 that window functions could systematically bias the estimated 
flux density of sources that are not at the tracking centre. 



value. In AMI, when the fringe rate drops below Vf ~ lO""* Hz, 
it becomes difficult to separate the signal from slowly-varying 
non-astronomical noise, so these samples are usually rejected. 
Some additional data will have to be rejected at low fringe rates 
because they cannot completely fill the PC data block. But the 
proportion of lost data is small for a telescope like AMI. 

4. Simulations 

We tested the oversampling concept with simulations of a model 
telescope based on AMI. The observing frequency is 15 GHz 
with a 6 GHz bandwidth, mixed down with an LO of 24 GHz 
to give an IF band of 6-12 GHz. Path compensation is applied 
in the IF at discrete steps. The shortest delay line is the same 
as the correlator's delay steps so that the gaps in the cross- 
correlation function can be filled over the course of a single path 
step. Fringe-rotation is applied in software, rather than in hard- 
ware so the fringes need to be sampled sufficiently fast. For a 
relatively short 5 m east-west baseline, we expect fringe periods 
longer than 50 s, so sampling at 1 Hz is more than sufficient. 8 
spectral channels are synthesised from the full 6 GHz bandwidth. 
We investigated both the real and complex correlators. The real 
correlator has 16 delay steps and the complex correlator makes a 
pair of measurements at each of its 8 delay steps. The parameters 
of the model telescope are summarised in Table [U 

Table 1. Model telescope parameters 



RF 

Bandwidth 

Synthesised sub-bands 

LO (Lower sideband reception) 

IF 

Dish diameter 

Baseline range 

Simulated baseline 

Integration 

Spectral channels 

Latitude 

Test source declination 



15 GHz (12-18 GHz) 

6 GHz 

0.75 GHz 

24 GHz 

6-12 GHz 

3.7 m 

5-20 m 

5 m east-west (250/1) 

Is 

8 (0.75 GHz each) 

52° 

52° 



Real correlator: 
Correlator delay steps 
Correlator delay steps 
Path compensator steps 
Tracking hour angle range 



16 samples 

83 ps (25 mm) 
83 ps (25 mm) 

±56 s 



Complex correlator: 
Correlator delay steps 
Correlator delay steps 
Path compensator steps 
Tracking hour angle range 



i complex samples 
167 ps (50 mm) 
167ps(50mm) 

±112s 



4.1. Simulations method 

We simulate the fringes measured at each delay step using 
Eq. (HI and (|7]i. We modelled an unresolved source at the centre 
of the field but the method will work for any sources because 
we are simply measuring the cross-correlation function. The an- 
tennas track the point source about the transit (when the source 
is due south) over an hour angle range of ±56 s for the real cor- 
relator or ±112 s for the complex correlator. This fills the CCF 
function over a range of At = +0.61 ns. The time-stream data is 
then fringe-rotated as outlined in Sect. 13. II The cross-correlation 
function is now a function of At and we grid the data points to 
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256 pixels (see the upper set of plots in Fig. |6]). We now have 
an oversampled measurement of the cross-correlation function. 
If the pixels are too narrow, some pixels may not contain any 
data points and this will degrade the spectrum. At the same time, 
the pixels must be fine enough so that the cross-coiTelation func- 
tion is not washed out. Individual pixels could be weighted by 
factors such as the number of data points that contribute to the 
pixel and the system temperature associated with the data points. 
Finally, we apply the fast Fourier transform (FFT) to the gridded 
data and extract the signal spectrum over the 6-12 GHz pass- 
band (lower plots in Fig. |6]l. The amplitude and phase can be 
caUbrated against an astronomical source by standard methods. 

4.2. Results and discussions 

4.2.1. Ideal correlator 

The recovered spectrum is in steps of 0.75 GHz and is the same 
as when the conventional one-shot method is used. The over- 
sampled spectrum extends from v,f = DC to 104 GHz and only 
a portion is shown in Fig.|6l The spectrum for the complex cor- 
relator in Fig. [6^ does not appear to suffer from aliasing. The 
arguments in Sect.[T] together with Fig.|2]suggests that the over- 
sampling method should reduce aliasing. 

In the real correlator, fringe-rotation introduces aliasing via 
the Fourier transform. This can be seen as small cycles in the 
data points (Fig. |6] and details in Fig. |2l). These arise from the 
time-domain alias cycles in Fig. [3] Fringe-rotation aliasing is 
worse near the edge of the cross-correlation function. The root 
mean square (RMS) deviation of the recovered spectrum from 
the ideal case is negligible (below 0.1 percent). But if there 
are delay errors, the RMS deviation may become significant 
and this will be discussed in Sect. 14.2.21 If using the oversam- 
pling method on real correlators, we would recommend hard- 
ware fringe-rotation. If software fringe-rotation is chosen, as 
shown here, the effects of fringe-rotation aliasing needs to be in- 
vestigated further. Gridding is equivalent to down-sampling and 
the conventional practice of introducing a lowpass anti-aliasing 
filter before gridding may be beneficial. With larger fractional 
bandwidths (Av/vrf), the small cycles from fringe-rotation alias- 
ing will become longer relative to the pixel size. 

4.2.2. Correlators with delay errors 

The oversampling method also corrects for delay errors preva- 
lent in analogue Fourier transform correlators. It is possible to 
determine the delay errors either by bench measurements or 
by tracking a source with t he path compensator held fixed (see 
[Harris & Zmuidzinasll2001l : H07). The delay errors are incorpo- 
rated into calculating At and gridded uniformly so the recovered 
frequency channels are also uniformly spaced. 

Figures [8^ and [8j) show simulations for both the real and 
complex correlators with delay errors of up to 10 percent of the 
delay steps. For the complex correlator, we assumed that the de- 
lay errors for each pair of detectors are the same. This is a rea- 
sonable assumption for our correlator architecture because we 
found that most delay errors arise in the long delay lines com- 
mon to both the real and imaginary signals (see H07). In both 
correlators, there are gaps in the samples from the irregular de- 
lay steps. The recovered spectrum is given by the spectral re- 
sponse of this distorted window function convolved with the true 
spectrum. This will inevitably worsen the inter-channel spectral 
leakage. 



The presence of delay errors makes aliasing resulting from 
fringe-rotation of the real correlator data noticeably worse 
(Fig.[8h and details in Fig.|7]l. This is consistent with our previ- 
ous findings (H07). When fringe-rotating the spectral channels 
from the real correlator, we assumed regularly-spaced centre- 
frequencies at the design frequencies. These may not be the best 
choice for an irregularly-spaced correlator and may even intro- 
duce small phase errors during fringe-rotation. Fringe-rotation 
aliasing gives rise to a low-level spread-spectrum noise in the 
amplitude plot (Fig. [8^). Perhaps the most significant evidence 
against applying software fringe-rotation to a real correlator is 
the potential for systematic errors. We found about 3 percent 
RMS deviation from the ideal spectrum when there are delay 
errors. The alias cycles seem to be largely responsible for this 
because the RMS deviation for a complex correlator with lag 
errors was less than 0.01 percent. 

It is worth noting that a correlator with delay errors will fun- 
damentally reduce sensitivity and no amount of signal process- 
ing can compensate for it. This is because the SNR of the recov- 
ered spectrum is optimised when sampled at Nyquist steps (or at 
half-Nyquist steps for complex correlators). Delay errors cause 
the noise components between the delay steps to be correlated, 
so the measurements are no longer independent and the SNR 
is degraded. The benefits of applying the oversampling method 
in the presence of delay errors are: (1) The recovered spectral 
channels are regularly gridded and are at the designed centre fre- 
quencies. (2) Normally, delay errors worsen the effects of alias- 
ing (H07). Fig.[8j5 and the RMS figure for the complex correlator 
suggest that the oversampling method reduces aUasing, even in 
the presence of delay errors. 



4.2.3. In the presence of noise 

In Figs.[8}; and[8]i, we added independent gaussian random noise 
to the fringes of a complex correlator. Noise in each time-frame 
are independent and it can be shown that in our idealised sys- 
tem, the noise between delay steps are also independent. The 
recovered spectra are noisier as expected but the oversampling 
method is robust in the presence of noise. It is also a linear pro- 
cess so successive data sets can be stacked (Fig. [8}l). Because 
the noise is independent between samples and also between the 
delay steps, the noise energy is spread over the whole of the 
oversampled spectrum. We only select spectral channels over the 
passband and discard the rest. This does not improve the SNR 
beyond what can be achieved by standard methods. The SNR 
improvement from applying this matched filter is the same as 
from averaging successive samples in the standard methods. 



4.3. Issues in a practical system 

Several issues will need to be addressed when using the over- 
sampling method in a practical system and design tolerances 
must be set. Firstly, more detailed simulations of the method 
should be conducted to quantify the reduction in aliasing. So 
far, we have assumed that the passbands of the detectors at each 
delay steps are identical. Any variation will restrict the effec- 
tive bandwidth and degrade the sensitivity. If there are variations 
in the passbands, each detector will trace out a different cross- 
correlation function because the cross-correlation function is re- 
lated to the Fourier transform of the passband. When the data 
points are merged to reconstruct the full cross-correlation func- 
tion, there will be discontinuities between the detector outputs. 
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However the degradation is probably comparable to standard 
methods. 

If the path compensators do not step exactly in multiples of 
one delay step (due to manufacturing tolerances in the path com- 
pensators), there will be unsampled gaps in At. Optimising the 
path compensation for an array can be more complicated and 
some delay ranges (At) may not be fully sampled. Alternatively, 
some samples may have to be rejected because they are cor- 
rupted. These missing samples will degrade the overall sensi- 
tivity in the same way as delay errors. Making the grid coarser 
may mitigate some of this but at the expense of washing out the 
cross-correlation. The data needs to be stable over the PC data 
block period. This could be worked round by opting for real cor- 
relators (with hardware fringe-rotation), which cuts down Tp, by 
a half. 
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4.4. Application to AMI 

Although both the real and complex correlators were trialled for 
AMI, the real correlator had been chosen (see H07). As demon- 
strated here, the oversampling method is not suited for real corre- 
lators relying on software fringe-rotation. Later on, it was found 
that the bottom two frequency channels in AMI (12-13.5 GHz) 
suffered from geostationary satellite interference. The effects be- 
came noticeable when observing at low declinations and it was 
decided to block this frequency range all together with hardware 
filters. Although this reduces the bandwidth of the instrument, it 
also goes some way towards reducing the aliasing on one side of 
the spectrum. This change avoided the need for the more com- 
plicated oversampling method and subsequently, we did not im- 
plement it in the AMI data reduction pipeline. 

5. Conclusions 

We have described a software-based method for overcoming 
aliasing in critically-sampled Fourier transform correlators. The 
method reduces aliasing by reconstructing the oversampled 
cross-correlation function from successive samples. Data from 
complex correlators can be fringe-rotated directly and efficiently 
in software. If using this method on real correlators, the signal 
should be fringe-rotated in hardware to avoid systematic errors. 
The edge spectral channels are usually discarded because 
they are the worst affected by aliasing. For an 8 -channel complex 
correlator simulated here, the sensitivity could be improved by 
up to 15 percent by retaining the edge channels. The additional 
computational overhead and complexity are offset by an im- 
provement in sensitivity and reduction in systematic errors from 
aliasing. This could be a significant benefit for Fourier transform 
correlators with less than 10 spectral channels. Some data with 
low fringe rates will have to be rejected, but for AMI-like tele- 
scopes with fractional bandwidths Av/vrf Z 0.4 and maximum 
baseline to dish ratios of D/d < 10, the fraction of data lost is 
relatively small. We have also shown that we can naturally com- 
pensate for delay errors to recover regularly-gridded frequency 
channels at the design frequencies. 
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Fig. 7. Details of the reconstructed cross-correlation function for 
the real correlator in Fig. |6h towards the edge of the delay range. 
The thick grey curve is the expected function from Eq. ^ with 
Tg set to 0. The dots are the data points generated from Eq. (fT2l l 
after the modified fringe rotation. The upper plot is for a perfect 
real correlator (from Fig. |6^) and the lower plot is with delay 
errors (from Fig. [8^). The solid vertical lines indicate the posi- 
tion of each delay in At at the beginning of the PC data block. 
The dashed vertical lines are the position at the end of the PC 
data block. Delay errors lead to missing samples or overlapping 
samples inside the narrow gaps indicated by these lines. Fringe- 
rotation aliasing is apparent in both cases and it is also worse 
when there are delay errors. 
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(a) The baseband signal (0-12 GHz) under and critically-sampled. The spectra show signs of aliasing. The cross-correlation function 
was sampled at 16 and 32 delays. 
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Fig. 2. The effects of different oversampling rates on the recovered spectra. The left-most plots show the cross-correlation function 
from Eq. (|7]i as a function of the residual delay At. This is a snapshot when the source transits an east-west baseline so the 
geometrical delay t„ = (case 1). The next set of plots to the right (case 2) show the cross-correlation function a little later when 
the centre of the envelope has drifted by a quarter of a delay spacing. The measurements at transit (case 1) are indicated by the 
grey circles and the measurements later on (case 2) are indicated by the black circles. The horizontal axis is the delay in At, which 
will be defined in Sect. 12. II The right-most plots are the amplitudes calculated from the delay measurements. The filled grey marks 
are at transit (case 1) and the dark circles are a while later (case 2). The two sets of marks should match but it can be seen in a) 
that they do not. This mismatch shows up as the temporal modulation seen in Fig. [3] In the upper row of a), although the 6 GHz 
signal bandwidth is critically-sampled at 16 delays, the 0-12 GHz baseband signal is undersampled. The two halves of the spectrum 
are side-by-side. The spectrum is aliased at both the high and low-frequency ends. In the lower row of a), the highest frequency 
component (12 GHz) is critically-sampled with 32 delays. Now only the upper frequency range suffers from aliasing. In b), the 
baseband signal is oversampled by a factor of 2 and 4 (sampled at 64 and 128 delays respectively). Now the two sets of spectral 
amplitudes match. Oversampling ensures a buffer between the two halves of the signal spectrum as well as with the image spectra. 
This reduces the spectral overlap and suppresses aliasing. In this case, oversampling the baseband frequency (0-12 GHz) by a factor 
of 2 is sufficient (an equivalent sampling frequency of 48 GHz at 64 delays). Oversampling above this does not reduce aliasing any 
further. As a minor technical detail, the forms of the spectra in b) are slightly different from the one in Fig. [T] Here (and also in 
Fig.O, we shifted the spectra by half a sub-band using the shift theorem of Fourier transforms. This ensures that the sub-bands span 
the desired frequency range and is discussed further in lHoIIer et al.l (l2007h . Subsequent simulations will assume channel-shifting. 
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Fig. 6. The oversampled spectral synthesis for a real correlator a) and complex correlator b). The top plots are the real components 
of the fringes. (1) The thick grey lines mark the ideal cross-correlation function we would expect from Eq. (01) for a) and the real 
component of Eq. ^ for b). The geometrical delay Tg for the thick grey curves are set to 0. The dark dots (seen as continuous 
curves at this resolution) are the modified-fringe-rotated data points measured by the correlator at each time sample (from Eq. ( fT2b 
and Eq. (fTTT i respectively). In b), the imaginary components have not been shown for clarity. The vertical lines indicate the range 
in At spanned by each delay; the curve in each section is traced out from the left boundary to the right. (2) In the second plot, the 
modified-fringe-rotated data points are pixelised into 256 pixels. For the complex correlator, the dark curve is the real component of 
Eq. ( fTTT i and the grey curve is the imaginary component. (3) The pixelised data are fast Fourier transformed. The bottom two plots 
show the amplitude and phase of the spectrum. The estimated spectrum extends to higher frequencies so only a portion is shown. 
The filled circles mark the synthesised spectral channels of interest over the physical passband of the system. 
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Fig. 8. a) and b): The oversampling method in the presence of delay eiTors. In both cases, the errors are up to 10 percent of the delay 
steps and there are visible gaps in the data. The solid vertical lines mark the left-most range of each delay and the dashed vertical 
lines mark the right-most range. The delay errors can be clearly seen from the relative alignment of these lines. The other effect in 
the real correlator a) is the higher level of fringe-rotation aliasing in the cross-correlation function (see details in Fig.|7]). The thick 
grey curves are the ideal cross-correlation curves from Eq. (|4]i with Tg set to 0. The dark points are the modified-fringe-rotated data 
points generated from Eq ( fTSl i and Eq. (fTTT i. c) and d): Gaussian random noise was added to the fringes of a complex correlator The 
instantaneous SNR is quantified by the ratio of the peak of the cross-coiTelation function to the standard deviation of the noise. In 
d), 16 data sets are stacked so that the recovered spectrum has the same SNR as c). The phases have been unwrapped arbitrarily 



